Fidelity decay in interacting two-level boson systems: Freezing and revivals 
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We study the fidelity decay in the fc-body embedded ensembles of random matrices for bosons 
distributed in two single-particle states, considering the reference or unperturbed Hamiltonian as 
the one-body terms and the diagonal part of the fc-body embedded ensemble of random matrices, 
and the perturbation as the residual off-diagonal part of the interaction. We calculate the ensemble- 
averaged fidelity with respect to an initial random state within linear response theory to second order 
on the perturbation strength, and demonstrate that it displays the freeze of the fidelity. During 
the freeze, the average fidelity exhibits periodic revivals at integer values of the Heisenberg time 
tn- By selecting specific fc-body terms of the residual interaction, we find that the periodicity of 
the revivals during the freeze of fidelity is an integer fraction of tn, thus relating the period of the 
revivals with the range of the interaction k of the perturbing terms. Numerical calculations confirm 
the analytical results. 
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I. INTRODUCTION 



Fidelity, or Loschmidt echo, was introduced to study 
the effect of perturbations on the dynamics of quantum 
systems, mainly in the context of quantum information, 
but also in the study of the stability of quantum sys- 
tems; for a review see [J]. Prosen and Znidaric noted 
that a freeze of fidelity occurs, i.e. fidelity will remain 
stable for very long times in the scale of the Heisenberg 
time, if the diagonal part of the perturbation matrix van- 
ishes 0. This lead to an additional interesting view on 
fidelity decay. In any given system, we can view any part 
of the Hamiltonian as the unperturbed system and the 
rest as perturbation. In particular in a many-body sys- 
tem the mean-field theory can be considered to be the 
unperturbed system, and the residual interaction as the 
perturbation. Fidelity thus serves as a measure of the 
quality of such mean-field approach. This was actually 
considered in Ref. Q in the context of a random matrix 
model using a two-body random ensemble, i.e. an en- 
semble that takes into account the two-body character of 
the interactions as well as the fermionic character of the 
particles. The result was rather surprising in the sense, 
that the freeze does not occur for the ensemble-average 
of fidelity, but it does occur for the median fidelity, or 
equivalently for the average of the logarithm of fidelity, 
also known as distortion [J, Q . In the bosonic case we 
can ask a similar question, and we have the pleasant sit- 
uation that such systems arc Liouville integrable in the 
semiclassical limit if the bosons are restricted to two lev- 
els This special case is attractive since it is experi- 
mentally accessible [rl-fioj. 



In the present paper we shall analyze this case in the 
framework of the embedded random matrix ensembles 
for bosons, i.e., we will allow in principle also three-body 
and higher-order interactions. Such ensembles have a 
long history for fcrmion systems [ll| - [l3j and a somewhat 
shorter one for bosons [14| - [l6| . More recently such en- 
sembles have even been defined for distinguishable parti- 
cles . For recent reviews see [l^ [l^ . For the bosonic 
case the freeze of fidelity is readily seen [23| , if we use the 
simplistic version of a mean-field theory that in addition 
includes the diagonal part of the residual interaction (in 
the representation in which the one-body component of 
the Hamiltonian is diagonal) . Besides the freeze of the fi- 
delity, we also uncover unexpected revivals at fractions of 
the Heisenberg time, which reflect the many-body resid- 
ual interactions. 

The paper is organized as follows: In Section [H] we 
define the fc-body two-level bosonic ensemble of random 
matrices and fidelity, for which the reference and per- 
turbed Hamiltonians are fully specified. In Section IIIII 
we calculate within the linear response theory the en- 
semble average of the fidelity and find that it is a Fourier 
cosine-series whose basic periodicity is the Heisenberg 
time. In Section |TV] we carefully select the perturbing 
off-diagonal fc-body terms and obtain fractional periods 
of the revivals during freeze, in units of the Heisenberg 
time. In Section |V] we summarize our work and outline 
the conclusions. 



II. DEFINITIONS 



The fc-body two-level bosonic ensemble of 
random matrices 



*On leave at Dept. de Matematica Aplicada i Analisi, Universitat 
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To define the fc-body Bosonic Embedded Ensemble 
of Random Matrices for bosons [l^ [13 1 consider n 



2 



spin-less bosons distributed over I single-particle states. 
These bosons are associated with the creation and anni- 
hilation bosonic operators and dj, respectively, with 
j = 1, . . . ,1. From here on, we shall focus on the two- level 
case 1 = 2, which is the simplest one, has some remark- 
able properties @, [U , and is also interesting from an 
experimental point of view for the two-component Bose- 
Einstein condensates [o, 10]. 

We denote the normalized two-level n-boson states 
as l/i) = \fi,n-n) = {M^^r^)-Ha\Yi4T~'"\^), where 
■^fij.,n = [mK*^ ~ m)']^^^ is a normalization constant and 
0) is the vacuum state. This is the occupation- number 
basis spanned hy fi = 0, ... n] the Hilbert-space dimen- 
sion is thus TV = n + 1. These states are coupled through 
a random fc-body interaction Hk{l3), which is written 

as illii 



k 

r.s—0 



V, 



Mr,k Ms,k 



(1) 



Here, k denotes the rank of the interaction, 1 < fc < n, 
and Vr^} are the fc-body matrix elements, which are inde- 
pendent Gaussian-distributed random numbers with zero 
mean and constant (fixed) variance Wq = 1. As in the case 
of the canonical random matrix ensembles [23| . Dyson's 
parameter /? distinguishes the cases according to time- 
reversal invariance: /3 = 1 is the time-reversal symmet- 
ric case, and /3 = 2 is the case where this symmetry is 
broken. Hence, the fc-body interaction matrix w'^^ is a 
member of the Gaussian orthogonal ensemble (GOE) for 
/3 = 1, or Gaussian unitary ensemble (GUE) for /3 = 2. 
Notice that Hk{(3) commutes with the number operator 
n = d|di -|- 0,202, i.e., the interaction conserves the total 
boson number n. 

As mentioned above, this ensemble for I = 2 presents 
some noteworthy properties: It exhibits non-ergodic level 
statistics in the dense limit [l6j . i.e., spectral or ensemble 
unfolding do not yield the same results when fc is fixed 
and n — > oo. In addition to this, for /3 = 1 the ensem- 
ble displays a large and robust quasi-degenerate portion 
of the spectrum for a wide range of fc, the Shnirelman 
doublets [i^l, while for /3 = 2 only seldom accidental 
quasi-degeneracies are observed These results are 

consistent with the fact that each member of the ensem- 
ble is Liouville integrable in the semiclassical limit Q. 

For later purposes, we write the fc-body Hamiltonian 
as HkW) = + VkiP), which are defined by 



Therefore, in the occupation-number basis, Hq^. (/?) is the 
diagonal part of Hk{l3), and Vfc(/3) contains only the off- 
diagonal contributions. 



B. Reference and residual interactions 

Fidelity is a measure of stability for small changes of a 
reference Hamiltonian [l[. Therefore, we must begin by 
defining the reference or unperturbed Hamiltonian Tig, 
and the perturbed one, which we write as Hx ~ Hq + 
XV, where A is the perturbation strength. We shall be 
interested in the case where the total interaction consists 
of a diagonal one-body interaction coupled to the fc-body 
two-level bosonic embedded ensemble. In particular we 
shall study the cases fc = 2 or fc = 3 which are physically 
the most relevant. 

To this end, we shall focus on a specific choice of 
and V, where we consider that V has zero diagonal ele- 
ments in the eigenbasis of the unperturbed Hamiltonian 
Ho; in this case, we shall refer to V as the residual inter- 
action. Notice that this case mimics the typical set-up 
in mean-field calculations, though it is also encountered 
in other cases of physical interest, e.g., when the per- 
turbation is a time-reversal symmetry breaking interac- 
tion [2^. We are interested in this type of reference and 
the residual interactions since they fulfill the conditions 
to observe the fidelity freeze , which implies longer sta- 
bility times. Therefore, we include the diagonal part of 
the fc-body interaction in the definition of the reference 
Hamiltonian, which we write as 

noW, A) = -^i^o.., + ^Ho, (/3). (3) 

In Eq. (|3]), we have normalized each term with the width 
of the spectrum Wk of the corres pon ding fc-body embed- 
ded ensemble, which is given by [l6l [l9l | 



W'^ = ^tT[Hkil3W ^ ^fik)+^-j^j2^Bi^-k), (4) 



N 



=0 



where 



(5) 



These expressions apply to the two single-particle level 
case {I = 2): the label B stands for bosons, the over- 
line indicates ensemble average, and A^^' (fc*) is the s-th 
eigenvalue of the ensemble-averaged correlation matrix of 
the bosonic fc-body embedded ensemble. 

We observe that, according to Eq. ([3]), the reference 
Hamiltonian depends explicitly upon A. Using the usual 
creation and annihilation rules, the unperturbed energy 
spectrum can be explicitly calculated 



Wi 



Wk 



Here, a > £2 without of loss of generality, and the coeffi- 
^2) cients Cj^^r are identically zero if r>/i or fc — r>n — /z. 



and otherwise are given by 



n — II 
k — r 



(7) 
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Then, the residual interaction consists simply of the re- 
maining off-diagonal matrix elements of the /c-body inter- 
action properly normalized by Wk, i.e., — Vk{l3)/Wk- 



Born expansion contains higher-order contributions with 
respect to A, since ^(^{t) does depend on A. Below, we 
shall restrict to the contributions up to second-order in 
A. 



C. Fidelity and fidelity amplitude 

Fidelity compares the time evolution of a given initial 
state under a reference Hamiltonian, with the time evolu- 
tion of the same state under a slightly different Hamilto- 
nian [ij. We use the Heisenberg time tn as the time unit, 
i.e., t = t' /tH, where tn = 2'Kh/d, and d = (ei — e2)/Wi 
is the mean-level spacing of 'Ho(/3, A). The unitary time- 
evolution associated with the reference Hamiltonian is 
given by Ul^{t) = rexp[-i27r/d?io(/3, A)t], where f is the 
time-ordering operator. Note that {t) inherits the A 
dependence from HoiP, A); we drop it from the notation 
to make it simpler. We denote by U^{t) the propagator 
associated with the perturbed Hamiltonian 'H\{I3). 

Considering an arbitrary initial state |^o)j the fidelity 
amplitude is defined as 



ft3,x{i) = (*o|M^,aWI*o) = (v&o|Wo^(-t)Wf (t)|*o), 
whose square modulus is known as the fidelity 



(8) 



(9) 



In ([U we have introduced the echo operator Mp^x^t), 
which corresponds to the time-evolution propagator as- 
sociated with the time -dep endent Hamiltonian Vi{t) in 
the interaction picture [2a] . 

Following [il, [2^, we use the Born expansion of 
Mp^xit), which we truncate at the second order; this 
approximation is referred as the linear response theory. 
Then, 



{^i\Mfl{t)\v) = 1 -iwfcA / dtiF/(ii) 

L Jo 



dii / dtaVf (ti)F/(t2) 1^^), (10) 



where Vf (t) = {—t)Vk{P)UQ (t) also depends on A and 
/?, and Wk ~ 2tt/ (dWk)- We notice that this second-order 



III. AVERAGING OVER THE ENSEMBLE 

We turn now to the calculation of the ensemble average 
of the fidelity amplitude. We use the occupation-number 
basis, where the reference Hamiltonian is purely diag- 
onal and the perturbation has diagonal elements equal 
to zero. We thus write the normalized initial state as 
1*0) = E^.^mIm)- Then, 

Jo 



t rt 



•wlX^Y.^l'^- / / di2(/x|vf(ii)vf(t2)k). (11) 



To carry out the calculation, we note that the terms 
on the r.h.s. of (jlip can be factorized in diagonal and off- 
diagonal contributions of the fc-body interaction matrix 
that is. 



{n\v,^{m = E^^^At) {^^\Vkm^^) . (12) 



(/i|y/(ti)^/(t2)|^) = E^M,p(*i)^P^(*2) (13) 



X {^l\Vkmp){p\Vkm^). 



where £^^^^{t) = exp [i(2^/d)(i?0(A) - i?°(A))t]. The 
time-dependence as well as the additional dependence 
upon A are contained in the factors p{t), which shows 
that they are determined by the /c-body diagonal ele- 
ments included in the reference Hamiltonian Ho{l3, A). 

Using the fact that the fc-body matrix elements have 
zero mean implies that Eq. (|12p is identically zero. For 
Eq. p3)) , we exploit the fact that the matrix elements of 
v^^-* are Gaussian distributed variables with fixed vari- 
ance, i.e., vl'^Jvl'^^ = {Sr,uSs,t + Sp^iSr,tSs,u)- For /3 = 2 
we obtain 



{„\Vkmp){p\VkmW) = S,^^Y.^flGl^Up~r\p-s)^gl%,^, (14) 

r^s 

I^JhWM = e^^[^27:[{^^-p)h + {p-l,)h]-\{wk\fY^^Gfl-^^^)^^ + ^^f)- (1^) 

r 

I 

Here, we denoted by tj^^p = X^r^^s G^^}rG^p}s{p — r\p — s), and used the fact that w/£=i(ei — £2) = 27r, which follows 




FIG. 1: Comparison between the ensemble-averaged fidelity 
given by Eq. (|16p (green/dashed curve) and numerical simu- 
lations (continuous/blue curve) for (a) /3 = 2 and (b) /? = 1. 
In this case fc = 2, A = 10"*^, n = 1024, ei = 0.7618036 and 
£2 = 0.9299698. Notice that the agreement extends to times 
beyond W^tH- The horizontal lines display the value of the 
fidelity freeze according to Eq. (|17|l . 



from the definitions of d and wi given above. For IJ- = p 
we have Q^l^ji ~ 0, since the matrix elements involved are 
exclusively off-diagonal. For the case (3 = 1, Eq. ((Ti)) has 
additional contributions which are of order 0{1/N), and 
shall be neglected. 

We note that, due to the first Kronecker delta in 
Eq. ([Tl)) . the integrand of the second term in pT|) de- 
pends only upon the time difference ti — ^2 • In addition, 
the second factor of Eq. (fTS)) is a correction due only to 
the diagonal fc-body interactions included in the reference 
Hamiltonian. Clearly, it is Gaussian in A, and therefore 
it contributes to corrections of the order of 0{X^). Con- 
sequently, this factor will induce corrections in the linear- 
response formula for the fidelity amplitude at least of the 
order ©(A"). 

We consider a normalized random initial state |5'o), 
and make the simplification ~ We carry 

up the time integrals and compute the square modu- 
lus retaining terms up to second order in A. Then, the 
ensemble-averaged fidelity to second order in A is given 




10"- 10^^ 10" 10' 10- 10-' 10-* IQ-^ lO'' 
t 

FIG. 2: Ensemble- averaged fidelity decay for various values 
of (a) the perturbation strength A, and (b) the number of 
particles n. Note that the time at which the fidelity freeze 
ends and decay sets back, t^, is independent of n. 



by 



4>) = ^- 



1 — cos[27r(^ — p)t] 



-0(A4). 
(16) 

Equation (|16p is the main result in this paper. In 
Fig. [l]we show the comparison of the ensemble-averaged 
fidelity predicted by Eq. (|16p with numerical calculations, 
both for /3 = 1 and /3 — 2. The results show excel- 
lent agreement even up to rather large values of t, when 
the fourth-order contributions in A eventually dominate 
and destroy the freeze of the fidelity. Notice that the 
agreement also holds for /3 = 1, in spite of the fact that 
Eq. (|16|) was obtained for /3 = 2. This confirms a pos- 
teriori that for (3 = 1 the corrections to ([H)) are a fac- 
tor l/{n+ 1) smaller as assumed, and can be neglected 
in leading order in the boson number. The fact that 
Eq. pop is not a perturbative expansion in time explains 
the good agreement for large values of t, which holds as A 
is small enough. In addition, for small times the results 
confirm the usual quadratic decay of the fidelity, namely, 

1 - i^Fl(t) (xt'^ + 0{t^). Fig. [2 displays the dependence 



of the ensemble-averaged fidelity with respect to A or n. 

More interesting and far reaching is the fact that fi- 
delity, up to second order in A is a Fourier cosine-series in 
t. The basic periodicity is precisely the Heisenberg time 
tn = 1, since the minimum difference in the occupation 
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numbers for the states |^) and \p) in one of the single- 
particle states is precisely 1. Indeed, the double sum over 
the basis states excludes the case /i = p, since for /i = p 
we have set ^^^^ = before the time-integration is car- 
ried up. This emphasizes the fact that the Fourier coefH- 
cients in Eq. (fT6|) are related to the off-diagonal residual 
interaction. Moreover, we observe that, in units of the 
iff, at integer values of time the time-dependent term 
of ([T5)) vanishes identically, and therefore revival of |^'o) 
are observed, as illustrated in Fig. [T] These revivals are 
not a full recovery of j^'p) though, since there may be 
corrections of higher order in A, and Eq. ([T6| has been 
averaged over the ensemble. We emphasize that the peri- 
odicity Tf = tn of the revivals of the ensemble-averaged 
fidelity for the bosonic embedded ensembles follows from 
the fact that ei ^ £2, which is responsible for the complex 
exponential factor of Eq. ([T5|) . 

Equation (jl6p permits to obtain an estimate of the 
freeze of the fidelity, which we denote by i^freeze- Con- 
sidering the minimum value for the fidelity when p6p 
is valid, that is, when the periodic revivals are observed 
because of the freeze of the fidelity, we have 



freeze 



= 1 - 



E 



(17) 



Equation (fT7|) predicts that i^frccze scales as A^. The 
scaling with respect to the number of particles is more 
involved since Wk and the coefficients Q^l^^p depend on 
n. Using Stirling's formula we obtain w\ ~ ri^^^'^ 



and 



^2k 



which yield the scaling 



where we took 



into account that the sum over the many-body states can- 
cels the normalization factor of the random initial state 
'I'o). These scaling laws are confirmed numerically as 
illustrated in Fig. |31 At this point we note that the time 
tf, during which the freeze of the fidelity lasts scales as 
tf. ^ A^ and is essentially independent of n; cf. Figs. [5] 
and[3l 



PERIODIC FRACTIONAL REVIVALS AND 
fc-BODY INTERACTIONS 



IV. 



As discussed above, Eq. predicts periodic time- 
revivals of fixed period Tf = tn ~ 1, independently of 
the rank fc of the residual interaction. This is a conse- 
quence of the fact that the residual fc-body interaction 
so far considered; contains terms which involve moving 
1,2, . . . , fc particles, from one of the single-particle lev- 
els to the other. That is, the many-body states coupled 
through the perturbation (/3) may differ at least in 
the occupation of one particle, and at most in the occu- 
pation of k. These differences are precisely the factors 
/i — p that appear in the Fourier expansion in Eq. (|16p. 
which denote the difference of number of particles in a 
given single-particle level. Clearly, the minimum differ- 
ence fixes the periodicity of the Fourier cosine series. 




FIG. 3: (a) Dependence of the freeze of the fidelity (circles) 
upon A, which confirms that 1 — Ffrcczc c>c A^. The right 
vertical scale indicates the time of duration of the freeze of 
the fidelity (data in squares) which shows the scaling tc oc 
1/71. (b) Dependence of -Ffreczc on n confirming the scaling 

1 ^freeze ^ ^ . 




FIG. 4: Ensemble-averaged fidelity decay for the residual 
interaction A^fc in Eq. (|18|) for k = 2 (green/dotted curve) 
and A; = 3 (blue/continuous curve), illustrating the fractional 
periodic revivals. These correspond to Tf — 1/2 and Tf = 
1/3, respectively. The magenta/dashed curve displays the 
result when all off-diagonal elements for = 3 of ([T]) are 
considered. 
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This explanation opens the following interesting pos- 
sibility. By selecting the actual perturbing terms within 
the A;-body residual interaction, we can actually tune the 
observed periodicity of the revivals during the fidelity 
freeze, in particular, making it different from Hcisenberg 
time Tf = 1. Indeed, we can fix the perturbation such 
that the only terms present move exactly k particles from 
one single-particle level to the other one. That is, we re- 
strict the general fc-body Hamiltonian ([1]) such that the 
off-diagonal terms are 

ICk - i^fc,o(aI)'X«2)'' + z;o,fc(4)''(«i)', (18) 

with Vk.o = (vo^k)* to ensure hermitecity. This choice of 
the residual interaction implies that the Fourier coeffi- 
cients vanish unless \^ — p\ = fc, that is, when the states 
l/z) and \p) precisely differ in the occupation of k parti- 
cles with respect to one mode. In this case, the argument 
of the cosine function in Eq. ()16|) is 27rfct, which implies 
that the periodicity of the revivals during the freeze of 
the fidelity becomes Tf = 1/fc in units of the Heisenberg 
time. Apart from the fact that this periodicity differs 
from 1, the important aspect is that the periodicity of 
the revivals during the freeze of the fidelity provides a 
direct measure of the rank of the interaction k of the 
residual perturbation ICk- These arc the fractional peri- 
odic revivals. 

This prediction is illustrated in Fig. HJ where we plot 
the ensemble-averaged fidelity for an interaction of the 
form (fT8|) for fc = 2 and fc = 3, and for comparison a case 
including all off-diagonal contributions of Hk in Eq. ([T]). 
The results clearly show Tf = 1/2, Tf = 1/3 for the first 
two cases, reflecting the value of the corresponding k- 
body interactions, and = 1, as expected from Eq. ()16|) . 

V. CONCLUSIONS 

In this paper, we have studied analytically and numeri- 
cally the fidelity decay in the fc-body embedded ensemble 
of random matrices for bosons distributed in two single- 
particle levels. We defined fidelity in terms of a refer- 
ence Hamiltonian, which is assumed to be in diagonal 
form, and a perturbed Hamiltonian which in addition in- 



cludes a purely off-diagonal residual fc-body interaction. 
This situation mimics the typical set-up in mean-field 
calculations, but appears also in other interesting physi- 
cal cases, such as time-reversal symmetry-breaking. This 
set-up fulfills the conditions to observe the freeze of the 
fidelity [3, , thus allowing for longer control of the sys- 
tem. 

We calculated the ensemble-averaged fidelity within 
the linear response theory up to second-order in the per- 
turbation parameter, which is a Fourier series whose ba- 
sic periodicity is equal to the Hcisenberg time tn- The 
analytical predictions are in good correspondence with 
the direct numerical results, confirming the presence of 
the freeze for the ensemble-averaged fidelity as well as 
the relevant scalings with respect to the strength of the 
perturbation and the number of particles of the system. 
The oscillatory part of this Fourier series cancels at inte- 
ger times of tn, thus manifesting the periodicity of the 
revivals. Selecting the off-diagonal terms of the fc-body 
residual interaction, in order that the actual perturbation 
couples only many-body states differing exactly by k par- 
ticles in the occupation number of either single-particle 
level, we showed that the periodicity of the revivals be- 
comes 1/fc in units of the Heisenberg time. Therefore, the 
periodicity of the revivals of the ensemble-average fidelity 
during freeze may be used as a direct measure to detect 
the rank fc of the perturbing interaction. This aspect 
may be interesting in the context of current efforts that 
address effects related to three-body interactions [13, HI] 
which are responsible, for instance, for atomic losses in 
ultra-cold bosonic gases. 
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